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ABSTRACT 



High frequency soft reverberation lags have now been detected from stellar 
mass and super massive black holes. Their interpretation involves reflection of 
a hard source of photons onto an accretion disk, producing a delayed reflected 
emission, with a time lag consistent with the light travel time between the irra- 
diating source and the disk. Independently of the location of the clock, the kHz 
Quasi-Periodic Oscillation (QPO) emission is thought to arise from the neutron 
star boundary layer. Here, we search for the signature of reverberation of the 
kHz QPO emission, by measuring the soft lags and the lag energy spectrum of 
the lower kHz QPOs from 4U1608-522. Soft lags, ranging from ~ 15 to ~ 40/xs, 
between the 3-8 keV and 8-30 keV modulated emissions are detected between 
565 and 890 Hz. The soft lags are not constant with frequency and show a 
smooth decrease between 680 Hz and 890 Hz. The broad band X-ray spectrum 
is modeled as the sum of a disk and a thermal comptonized component, plus 
a broad iron line, expected from reflection. The spectral parameters follow a 
smooth relationship with the QPO frequency, in particular the fitted inner disk 
radius decreases steadily with frequency. Both the bump around the iron line in 
the lag energy spectrum, and the consistency between the lag changes and the 
inferred changes of the inner disk radius, from either spectral fitting or the QPO 
frequency, suggest that the soft lags may indeed involve reverberation of the hard 
pulsating QPO source on the disk. 

Subject headings: Accretion, accretion disks. X-rays: binaries, stars: neutron, galaxies: 
active, X-rays: galaxies 



- 3 - 



Introduction 



Kilohertz quasi-periodic brightness oscillations (kHz QPOs at frequ ency above ~ 300 



Bradt et al. 



2000 



19931) 



for a 



Hz) have been observed with the Rossi X-ray Timing Explorer ( RXTE, 
from about thirty neutron star low-mass X-ray binary systems (Ivan der Klis 
review). Those signals have triggered much excitements, because their frequencies match 
the orbital frequencies of matter orbiting very close to the neutron star, hence may be 
related to strong field general relativity effects. Most of the analysis performed so far on 
kHz QPOs has focussed on their frequency, RMS amplitude, quality factor, relation to 
continuum spectral parameters and variability at lower frequency, but very little is known 
about their energy dependent time lags. Soft lags were however repor ted using early RXT 



observations of the lower kHz QPOs of 4U1608-522 and 4U1636- 536 (IVaughan et al. 



Kaaret et al. 



1999 



Markwardt et al. 



20001 : 



de Avellar et al. 



1998 



201l|). In 4U1636-536, at - 830 



Hz, the soft emission in the 3.8- 6.4 keV band lags the hard emission above 9.3 keV by 



about 25 fis (iKaaret et al 



19991 ). a value that is consistent wi th the one measured from 



4U1608-522 at about the same frequency 
the 4-6 and 11-17 keV bands). Recently 



(Vauehan et al. 


1998 


) 


'~ 27fis at 830 Hz between 
) have extended on those 


de Avellar et al 


( 


2013 



previous analysis by a more thorough study of both sources. 

The detection of soft lags from kiloHz QPOs of neutron stars has to be placed 



in a broader context, where soft lags are now being routine - 



frequency variabi 



De 



Marco et al 



ity o f stellar mass black holes 



2011 



Emmanoulopoulos et al. 



Utt 



201 



ey et al 



V det ected in the high 



20111). activ e galactic nuc 



Cackett et a 



20 



Fabian et al 



ei 



20131 ). Those lags have 



20131 ) and even recently from cataclysmic variables (iScaringi et al. 
great significance as they probe the dynamics of the very inner regions of the accretion flows 
onto compact objects. They are a natural prediction of the reflection model, in which a hard 



X-ray source irradiates the accretion disk, producing a delayed emission with a reflection 



-4- 



spectrum (with prominent Iron Ka line eraission), that is skewed by relativistic effects 



close to the central object (see 



Miller 


2007; 


Reynolds 


2013 



for reviews). The magnitude 



of the lags is consistent with the light travel time between the irradiating source and the 



reflector. T 



le detection of a lag in the Kg line em is sion provides clinch ing evidence for such 



a scenario ( IZoghbi et al. 



2012 



Kara et al. 



2018a 



Zoghbi et al. 



20181 ) 



Independently of the mechanism setting up the frequency, it is generally agreed that 
the QPO modulated emission arises from the boundary layer between the neutron star 



surface and the accretion disk (iGilfanov et al 



2008 



Abramowicz et al. 



20071 ). The strongest 



evidence comes from Fourier frequency resolved spectroscopy showing that the QPO energy 
spectrum resembles the energy spectrum of the boundary layer. The QPO spectrum is 
significantly harder than the average continuum emission, and certainly the disk emission, 



also detected in the energy spectrum (iGilfanov et al 



20031) 



In this paper, we investigate wether the soft lags detected in kHz QPOs can be 
interpreted in the framework of the reverberation lags recently detected in AGNs and stellar 
mass black holes. For this purpose, we measure the lags and the spectral parameters of the 
continuum emission as a function of the lower kHz QPO frequency. We limit the present 
analysis to the highest quality RXTE data of 4U1608-522, showing strong kHz QPOs 
between 565 and 890 Hz, as to maxi mize the sensitivity to lag detections (basically the 



same data used in 



Barret et al 



de Avellar et al. 



20051). We note that our work extends on 
( l2013l ). who did not use data above 800 Hz for 4U1608-522 and did not combine the lag 
measurements with the spectral analysis of the continuum emission. In the next section, we 
describe the technique used and the main results. 
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2. Data analysis 

To measure the time lags, we have first extracted events in two adjacent energy bands: 
3-8 keV and 8-30 keV, using high time resolution science event data. The cross spectrum of 
the soft and hard b and time series gives the phase difference between the two energy bands 



( jNowak et al. 



19991 ) ■ This phase difference is converted back into a frequency-dependent 



time lag: r(z/) = A0/27rz/. In this paper, a positive lag means that the soft band lags behind 
the hard bancjl]. We have computed power density and cross spectra over segments of a 
fixed duration T and have averaged M of those to produce an individual Cross Spectrum 
(CPj). The lower kiloHz QPO frequency {vi^i) can vary by tens of Hz over the typical 
duration of the observation (~ 3000 seconds). The analysis presented below applies only to 
the lower kHz QPO of 4U1608-522, and is not suited for the upper kHz QPO, which is much 
harder to detect (mostly because of its larger width). As to maximize the lag detection 
sensitivity, we wish to bin each individual CPj over adjacent frequencies, sampling the QPO 
profile. This implies knowing the QPO centroid frequency and its width {wi^i, FWHM) in 
each CPj. In order to do this, we have first reconstructed the time evolution of the QPO 
frequency, extracting Fourier PDS using events in a 3-30 keV reference band. The reference 
energy band chosen maximizes the signal-to-noise ratio of the QPO detection, and therefore 
enables a better determination of the QPO frequency than in the soft and har d bands, vi^i 



20121 ). We have 



was obtained with a maximum likelihood technique (iBarret fc Vaughanl 
also shifted-and-added all the PDS to the mean QPO frequency to obtain the mean QPO 
width {w) over contiguous segments of data. Here we have assumed wi^i = w and averaged 
the CPj over a frequency bandwidth of twice w. As a final product of the analysis, for each 
observation, we have frequency dependent time lags integrated over segments of M x T 



^This is the opposite convention to the one that is commonly used in AGN studies, where 
soft lags are usually negative. 
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second duration and frequency bandwidth 2 w. Hereafter, we have assumed T = 4 sec, 
which provides adequate frequency resolution to sample to the QPO profile, whose width 
varies from ~ 3.5 to ~ 10 Hz in the data considered here. 

In the second segment of the March 3rd, 1996 observation (ObsID 10072-05-01-00), the 
lower kHz QPO can be significantly detected in the soft and hard X-ray bands on timescales 
as short as 32 seconds, as a combination of a large source count rate (the 5 detector units 
of the proportional counter array were operating) and high coherence. In Figure [H we show 
the soft lags between the 3-8 keV and 8-30 keV bands, measured on 32 and 128 seconds 
respectively (i.e. M = 8 and 32). As can be seen, increasing M yields fully consistent 
results, while the scatter on the lags and associated error bars have been reduced. Note 
that the lag is consistent with being constant in both cases, a round a mean yalue o f ~ 23/xs. 



Such a value is fully consistent with previous measurements (IVaughan et al. 



19981 ). 



2.1. Frequency dependent time lags 

Across its full frequency span, the QPO cannot be detected on timescale as short of 32 
seconds in the two adjacent bands. For this reason, we have adopted 128 second integration 
time for the lag measurements (i.e. averaging M = 32 intervals of 4 seconds). In Tabled] 
we list the observations considered in this paper, together with the mean lag measured. We 
have split the data in ObsIDs first, but also within an ObsIDs in the contiguous segments of 
data provided by the science event files. In deriving the mean soft lag, we have considered 
only segments of 128 seconds, in which the QPO is detected above ~ 4.5cr (excess power) 

R > 1.5, wher e R is 



Boutelier et al. 



fl2010h 



in both the soft and hard energy bands (this corresponds to a cut o 
the ratio between the Lorentzian normalization and its la error, see 
for a discussion). For a given integration time, the QPO significance depends on the source 
count rate and frequency; the latter parameter defining the QPO RMS amplitude and 
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width (IBarret et al. 



20051 ) ■ In the dataset considered, the significance of the QPO goes 
down at both ends of the frequency range, more dramatically at the lower end where the 
source count rate was the lowest. Nevertheless, adopting the cuts above, enables to detect 
the lags between ~ 565 and ~ 890 Hz. In Figure [21 we show the soft lags measured on 
128 seconds as a function of the QPO frequency. The lags have also been binned into ten 
adjacent frequency bins. A clear trend is now obvious in the data, in particular above 680 
Hz where the soft lags smoothly decrease with frequency. 



2.2. The energy spectrum of the lags 

We now wish to examine how the lags vary with energy. In each energy bin, using the 
same procedure as described above, the lag is computed between the light curve in that 
bin and the light curve in the reference energy band (3-30 keV), where the signal to noise 
ratio of the QPO detection is the highest. The light curve in the energy bin considered is 
subtracted from th e reference light curve, to ensure that Poisson noise remains uncorrelated 



flUttlev et al. 



20111 ). For each segment of 128 seconds, hence for each QPO frequency Ui, 
for each energy bin, we thus have a measurement of the lag, as the average of M = 32 
values. A mean lag-energy spectrum can be computed by averaging all the measurements 
weighted by their errors, within a continuous science event file duration, and even within 
one ObsID, provided that the channel energy boundaries do not change. Two examples of 
such lag spectra are shown in Figure in [3] for the first two ObsIDs which provide the better 
statistics (10072-05-01-00 and 30062-02-01-000). Note that the data were unfortunately 
not recorded in the same spectral mode in the two ObsIDs: 64, 32 channels. Although of 
limited significance, for the first ObsID where the source count rate is the highest, it is 
interesting to note the lag spectrum is fairly steep below 4 keV, and seems to fiatten around 
5-8 keV, where the weak and broad Iron line is significantly present in the energy spectrum 
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(see below), before steepening again. There is also a dip in the lag spectrum just around 
6.4 keV. This could be due to the presence of a constant and narrow iron line component 
(e.g. the core of the broad line), produced at larger distances from the neutron star or over 
a wider range of radii than its red wing. The same features are however not so clear in the 
other segment of data shown in Fig. [3] (right panel), although the data were recorded with 
a poorer spectral resolution, which could easily hide the structure seen at higher spectral 
resolution. 



2.3. Time averaged continuum emission spectrum 



In order to place the results described above into context, it is worth looking at the time 
averaged energy spectrum to see how it changes with respect to the Q PO frequency. One 



of the most extensive spectral analysis of the RXTE data, presented in 



Gierlihski fc Done 



( 120021 ) sampled the full set of spectral states of the source (it did not include any timing 
information though). When grouping data based on spectral colors (as to improve the 
statistics), they showed that across the Z-shaped color-color diagram, all X-ray spectra 
could be represented as the sum of a disk blackbody component, plus a comptonized 
component and a weak reflection component from an ionized disk. Both the reflection 
and its self consistently evaluated iron line emission were statistically detected in all the 
spectra. In this paper, we have extracted Standard 2 spectra for all the time intervals of 
Table [H We have combined the data from all PCUs units and used the latest calibration 
data available in heasoft 6.12, adding a systematic unc ertainly of 0.25% to the data. 



1989h . Although 



We have frozen the column density to 1.5 x 10^^ cm ^ (jPenninx et al. 
there is rnore than one possible decom position of the X-ray spectrum of 4U1608-522, (e.g. 



Lin et al. 



2007 



Takahashi et al 



Gierlihski fc Pond (|2002[ ) (see also, e.g 



201ll ). we have assumed the widely used spectral model o f 



Barret et al 



2000 



Barret 



2001 



Di Salvo fc Stella 



20021 : iTarana et al.l 120081 : iCackett et al.l 120101 . for a disc ussion of models). Howe ver, with 
the reduced statistical quality of our data compared to iGierlihski fc Pond ( 120021 ). we were 



unable to obtain stable fits when adding the self-evaluted reflection of the Comptonized 
emission (this is largely a consequence of having two curvy components and a broad iron 
line profile within the relativel y limited bandpa ss of the PGA). We have also used the 



Comp tt comptonization model (iTitarchuk 



19941 ) instead of the nthcomp model (IZycki et al. 



19991 ). We have checked however that both models provide consistent fits. The addition of 
a diskl ine component imp royed the fit significantly in all segments of data, as previously 



found (iGilfanov et al 



2003 



Gierlihski &: Done 



20021). Unfortunately, the data do not allow 



to constrain the inner disk radius of the diskline model, which we have frozen to lORg (for 
a disk inclination of 45 degrees). Different values (e.g. QRg) would fit the data as well. The 
normalization ( Ndrr) of the diskbb cora ponent was converted into an inner disk radius 
(Rin), following iGierhhski fc Donel (120021 ): 



Rj 



O.QIN 



1/2 
DBB 



2.7 



D 



fcol 



0.5 



cos I 



1/2 



km 



(1) 



7] J \3.6 kpc 

where D is the source di stance, i is the disk inchnat ion angle, fcoi is the ratio of the color 

"actor for the 



to effective temperature (IShimura &: Takahara 



19951 ) and t] is the correction 



inner torque- free boundary condition [rj = 2.7 for R^ = 6 -Rj) ( iGierlihski et al. 



19991 ). 



The best fit spectral parameters are listed in Table |2] and plotted in Figure H] against 
the lower kHz QPO frequency. As can be seen, all spectral parameters show a smooth 
behavior with frequency, suggesting that our spectral decomposition is robust. From table 
2, the significance of the line can be inferred from looking at the with and without the 
line. In our observations, the Ax^ varies between ~ 15 up to ~ 60 for 2 additional degrees 
of freedom (and a total of degrees of freedom around 40), making the diskline component, 
highly significant. 
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3. Summary of the results 



As advocated for instance by 



Uttley et al. 



(1201 ll ). frequency and energy dependent 



time lags give additional information, compared to time average energy and power density 
spectra, providing new and independent constraints on models for the emission mechanisms 
and the location at which they take place in the system. Our main results can be 
summarized as follows: 

• The soft lags between the 3-8 and 8-30 keV QPO photons vary with frequency, 
dropping at both end of the frequency range investigated (from 565 to 890 Hz). In 
particular, the lags decrease smoothly from about 40 fis to 15 fis, while the frequency 
varies from 680 Hz to 890 Hz. Our results are consistent with the values reported by 



de Avellar et al.i (|2013|) below 800 Hz. The same curvy shape in the l ag-frequency plot 



seems to be present in the data of 4U1636-536 (see Figure 3 top panel 



de Avellar et al. 



20131 ). suggesting that this may be a common feature of lags from lower kHz QPOs. 



It should be stressed however that this is not the intrinsic lags that is measured, 

y to contribute to the 



Zoghbi fc Fabian 



201lh . 



due to the fact primary direct and delayed emissions are like 
emission in both energy bands (the so-called lag dilution, e.g. 
With this caveat in mind, it is still worth trying to relate the change in the lags with 
frequency with proxies of the disk location, such as the kHz QPO frequencies or the 
inner disk radii inferred from the spectral fitting of the continuum emission. 

• When the lower kHz QPO frequency varies from ~ 600 to ~ 900 Hz, the upper kHz 
QPO frequency varies from ~ 900 to ~ 1070 Hz. For a neutron star of 1.4-2.0 Mq, 
assuming that either frequency is an orbital frequency at the inner edge of the disk, 
this implies a change of the inner disk radius of ~ 6.2 — 7.0 km (or ~ 3.5Rg for a 
1.4 Mq canonical neutron star) or 2.2 — 2.4 km, or equivalently a light travel time 
difference of 23-21 or 7-8 /xs. These values are not very different from the span of 



-li- 



the soft lags measured. Interestingly, considering only the decreasing part of the 
lag-frequency variation of Fig. [T] (between 680 and 900 Hz), the relative change of 
the lags by ~ 25yus would be very comparable to the light travel time difference, 
associated with the relative change of the orbital radius inferred from the lower kHz 
QPO frequency change (4.7 km or ~ 16 /xs). 

• The spectral analysis of the continuum emission is consistent with a picture in which 
the disk gets closer to the neutron star, while the QPO frequency increases. The 
relative change of the inner disk radius inferred from fitting the disk component is 
again broadly consistent with the span of the lags measured. It is also consistent 
with the relative radius change assuming that the lower kHz QPO is an orbital 



frequency. The distanc e to the sour ce is not very well 



( Nakamura et al. 



1989I ). as used by ( 



recent value of 5.81^° kpc (iGiiver et al 



Gierlihski fc Done 



mown, ranging from 3.6 kpc 



2OO2I ) in equation 1, to a more 



Suleimanov et al. 



20101 ) from type I burst analysis (see also 



2OIII )). Assuming for instance a distance of 4.4 kpc or 5.0 kpc in 



equation 1, for a neutron star of 1.4 Mq or 2.0 Mq respectively, would provide a 
perfect match between the inferred inner disk radius variations from spectral fitting 
and the orbital radius changes assuming the lower kHz QPO is providing the orbital 
frequency (see Fig. [5]). On the other hand, this would not work if one assumed that 
the upper kHz QPO is providing an orbital frequency, because the change in radii 
associated with the change of frequency would be too limited (~ 2 km as opposed to 



6 — 7 km). We note that the debate o n which one of the two kE 



is an orbital frequency is not yet settled (lOsherovich fc Titarchuk 



z QPO frequencies 



1999 



Sanna et al. 



2OI2I ). althoug h most mode 



freque ncy (e.g. 



Miller et al. 



s predict that it is 



1998 



van der Klis 



:he u pper QPO that pro vides an orbital 



20001 ). As discussed by 



Lamb fc Miller 



( 120011 ) ■ it is also possible that the observed upper kHz QPO frequency is significantly 
lower than an orbital frequency, which may then span a wider frequency range and 
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reach higher frequencies, so that the associated change of radius would be actually 
larger than observed. In absolute units (km), to make it consistent with the inner 
disk radius changes inferred from spectral fitting (Fig. H]), one would need to tune 
some of the parameters in equation 1, so that Rin drops by some 20% or so. This does 
not seem unfeasible given the uncertainties on the four parameters of equation 1 (e.g. 
the source distance). As a side note, interpreting the drop of the quality factor of the 
lower kHz QPO, assumed to be or bital, as a signatur e of the innermost stable circular 



orbit would then be more natural (IBarret et al. 



20061 ). However, it would be necessary 



to explain how the system can generate frequencies higher than the ISCO frequency 
(i.e. the upper kHz QPO). This would also imply large masses for the neutron star 
in kHz QPO sources. In the case of 4U1608-522, if the ISCO frequency is at 900 
Hz (the highest frequency detected here), then the gravitational neutron star mass 
would be around 2.6 Mq, a ssuming the dimen sionless angular momentum of the star 



(j = cJ/{GM^)) to be 0.1 flMiller et al. 



19981 ). Such a value would b e significantly 



large r than the mass of the most massive neutron star known to date (jPemorest et al. 



20101 ). 



• Evidence for reflection comes from the presence of a broad i onized Iron line, significan tly 



detected in all spectra, although weak, as pointed out by iGierlihski fc Pond (120021 ). 



This is a necessary ingredient for the reverberation scenario. The presence of several 
components (from the disk, the comptonization cloud and the reflection) in the PGA 
bandpass, together with its limited spectral resolution do not enable us to extract the 
parameters of the reflection component (e.g. ionization state of the reflector). There 
is a trend for the equivalent width of the line to decrease with frequency that needs 
to be understood. 



• Understanding the lag energy spectrum would require detailed simulations of the 
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transfer function of the irradiated disk, but looking at Figure [2l it is quite striking that 
they may aheady contain some very valuable information. This is particularly true in 
the highest count rate observations (recorded with the highest spectral resolution), 
which show both a soft excess and a broad bump around the broad iron line, present 
in the averaged energy X-ray spectrum. It is thus tempting to interpret the soft 
excess as thermal reprocessing of the hard X-ray photons onto the disk, and the bump 
around the broad iron line as the signature of delayed reflection. Further support 
for this interpretation comes from looking at the energy spectrum of the lag [21 The 
mean lag around 15 keV is -20 fis while the lag around 5-7 keV is about 10 fis. The 
total lag of 30 fis corresponds to a distance of ~ 4Rg for a 1.4 M© neutron star. For a 
neutron star of 10 km or 5Rg radius, this would put the inner disk radius at ~ 9Rg; 
a value that would fit the data with the diskline model and that is t ypically inferred 



from fitting the broad iron line of neutron stars (iCackett et al. 



l]) 



20101). More detailed 



modeling is still required to test the reverberation scenario. 



Conclusions 



To conclude, although alternative explanations for the soft lags have been proposed, 
for instance, as being due to a n intrinsic spectral softening of the emission along the 



QPO cycle (IKaaret et al. 



19991 ). or due to temperature oscillations in th e corona driving 



with some delays the 



de Avellar et al. 



temperature oscillations of the soft photon source (ILee et al. 



2001 



20131 ). the recent detection of soft lags in a variety of compact objects, 
from stellar mass black holes to AGNs, and the evidence accumulated above suggest 
that the lags measured involves reverberation of the hard pulsating source located in the 
boundary layer and interacting with the accretion disk. It is worth noting that for a similar 
geometry, down scattering of hard X-rays (from the accretion column) in a cool medium 
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(the accretion disk) was proposed t o explain the ~ 100/^5 soft lags observed in the pulsed 



emission of millisecond pulsars (e.g. 



Falanga fc Titarchuk 



20071 ). In our data, the distance 



between the irradiating hard source and the reflector varies in a way that can be tracked by 
the frequency of kHz QPOs. 

Our results illustrate clearly the power of lag measurements to probe the innermost 
regions of the accretion flows onto compact objects. Even if the surfaces and faster 
variability timescales of neutron stars make them more complex to analyze than black 
holes, kHz QPO sources offer stable and coherent clocks that allow us to measure time lags 
comparable to the light crossing time of a few Kg (the light crossing time at IRg is about 
7fis for a canonical neutron star), just as for AGN. Interpreting the detailed shape of the 
lag energy spectra, and determining intrinsic lag values will require modeling the transfer 
function of the response of the accretion disk to the irradiating source, which is beyond 
the scope of this paper. Putting a yardstick on a neutron star accretion disk will then be 
within reach. The joint timing and spectral analysis presented here should also be extended 
to other source s for which kHz QPOs were detected cleanly with the RXTE PGA (such 



as 4U1636-536, 



de Avellar et al. 



20131 ). The discovery of patterns similar to what we have 



found for 4U1608-522 (e.g. broad bump around the iron line in the energy spectrum of the 
lags) would add support to the reverberation scenario. It would also prepare the ground 



for breakthrough 



1 observations to be performed by the next generation of timing missions. 



such as LOFT (IFeroci et al. 



20111 ). whose combination of large effective area and improved 
spectral resolution will uncover exquisite details about the soft lags for both the lower and 
upper kHz QPOs, providing errors on the lags twenty time smaller and enabling truly iron 
line reverberation studies. This holds great potential to constrain the neutron star mass 
and to understand the underlying physics of the QPO emission and its propagation in the 



strong field region around the neutron star. 
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Fig. 1. — Soft lags (in fis) measured between the 3-8 keV and 8-30 keV bands for the 
second March 3rd, 1996 observation of 4U1608-522, for integration time of 32 (bottom) and 
128 seconds (top panel) respectively. In both cases, the soft lags are consistent with being 
constant with a mean value around ~ 23/is. The best fit value is represented by the dashed 
line. 
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Fig. 2. — Top) Soft lags between the 3-8 keV and 8-30 keV time series, binned in 10 adjacent 
QPO frequency intervals. The lags show a clear dependency with frequency, with the most 
significant trend being the smooth decrease at frequencies above ~ 650 Hz. Bottom) Same 
as above, but the individual soft lags measured on 128 second timescales are shown. 
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Fig. 3. — Mean soft lag energy spectra measured for the ObsIDs 10072-05-01-00 (left) and 
30062-02-01-000 (right). The lag-energy spectrum on the left is the best that could be 
obtained, as a combination of a narrow spectral binning and high count rate statistics. Note 
in particular the broad bump between 5 and 7 keV, suggestive of iron line reverberation. 
The dip around 6.4 keV may be associated with a constant and narrow iron line component 
(e.g. the core of the broad line), produced at larger distance from the neutron star or over a 
wider range of radii than the red wing of the line. It is also worth noting that the mean lag 
around 15 keV is -20 fis while the lag around 5-7 keV is about 10 fis. The total lag of 30 
fis corresponds to a distance of ~ 4Rg for a 1.4 Mq neutron star. For a neutron star of 10 
km or 5Rg radius, this would put the inner disk radius at ~ 9Rg] a value that would fit the 
data with the di skline model and th at is typically inferred from fitting the broad iron line 



of neutron stars (ICackett et al. 



20101 ). However, the trend seen on the left hand side plot is 
not so obvious on the plot on the right hand side, but note that the data were recorded with 
a much lower spectral resolution. 
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Fig. 4. — Variation of the best fit spectral parameters as a function of the lower kHz QPO 
frequency. Left) From top to bottom, the iron line equivalent width (in keV), the inner disk 
radius (in km), computed from equation 1 (with canonical parameters), and the inner disk 
temperature (in keV). Right) From top to bottom, the optical depth, the electron temper- 
ature (in keV) and the seed photon temperature (in keV). Errors on the fitted parameters 
are given at the la level. The error bar on the frequency is determined by the spread of 
QPO frequency over the continuous science event file. As can be seen, all best fit spectral 
parameters, but the seed photon temperature (at the highest frequencies), show a smooth 
behavior with frequency. 
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Fig. 5. — Same as Figure H] but the inner disk radius inferred from spectral fitting is com- 
puted from equation 1 for a distance of 5 kpc. The sohd hne represents the orbital radius 
derived from the orbital frequency (assumed to be provided by the lower kHz QPO) around 
a neutron star mass of 2 Mq. The Spearman correlation coefficient is -0.69, corresponding 
to a null hypothesis p robability of ~ 5 x 10^^. The source distance is 5.8lfg kpc according 



to 



GuveretaL 



(120 loh . 



Table 1. Log of lags detections for 4U1608-522 on timescale of 128 seconds. The table lists the observation date, the 
ObsID, the number of the science event file within the ObsID, the start date of the observation data set, the 
frequency span, the mean source count rate, the mean time lag and associated error over the science event file 
considered, the and the number of degrees of freedom for a fit by a constant, and finally the exact energy bands 
used for the soft band (the first two energies), the hard band (the second and third energies). 



ObsID n Start date Frequency space Count rate Mean lag (d.o.f) Energy bands 



10072- 


-05- 


-01- 


-00 


1 


1996/03/03-19:17:21 


847.2-892.5 


2943.3 


14.3±3.3 


21.5 


(24) 


3.0-8, 


.0- 


-30.5 


10072- 


-05- 


-01- 


-00 


2 


1996/03/03-20:53:19 


823.1-836.3 


2827.9 


22.6±2.4 


24.2 


(27) 


3.0-8, 


.0- 


-30.5 


10072- 


-05- 


-01- 


-00 


3 


1996/03/03-22:29:18 


829.3-840.4 


2861.9 


25.1±4.3 


19.3 


(13) 


3.0-8, 


,0- 


-30.5 


30062- 


-02- 


-01- 


-000 


1 


1998/03/24-15:36:34 


750.7-772.0 


1954.0 


27.7±3.3 


16.4 


(20) 


2.9-8, 


,2- 


-32.4 


30062- 


-02- 


-01- 


-000 


2 


1998/03/24-17:06:23 


796.4-819.1 


2042.9 


23.9it3.8 


36.1 


(26) 


2.9-8, 


,2- 


-32.4 


30062- 


-02- 


-01- 


-000 


3 


1998/03/24-18:42:22 


735.7-776.3 


1943.7 


25.2±3.5 


27.8 


(26) 


2.9-8, 


.2- 


-32.4 


30062- 


-02- 


-01- 


-000 


4 


1998/03/24-20:18:20 


665.5-709.3 


1818.4 


38.0±3.5 


19.1 


(25) 


2.9-8, 


.2- 


-32.4 


30062- 


-02- 


-01- 


-00 


1 


1998/03/24-22:01:31 


643.7-667.2 


1751.6 


26.1±4.9 


19.5 


(23) 


2.9-8, 


.2- 


-32.4 


30062- 


-02- 


-01- 


-00 


2 


1998/03/24-23:53:30 


643.2-665.6 


1415.0 


40.3±8.6 


19.4 


(15) 


2.9-8, 


,2- 


-32.4 


30062- 


-02- 


-01- 


-01 


1 


1998/03/25-13:54:37 


593.5-607.9 


1547.6 


26.7±9.4 


15.3 


(18) 


2.9-8, 


,2- 


-32.4 


30062- 


-02- 


-01- 


-01 


2 


1998/03/25-15:30:25 


565.9-580.0 


1531.1 


12.2±29.5 


12.1 


(8) 


2.9-8, 


,2- 


-32.4 


30062- 


-02- 


-01- 


-01 


3 


1998/03/25-17:46:60 


567.1-576.3 


1467.5 


ll.litll.l 


1.3 (6) 


2.9-8, 


,2- 


-32.4 


30062- 


-02- 


-01- 


-01 


4 


1998/03/25-18:42:22 


595.4-631.4 


1512.4 


29.3±8.9 


10.5 (13) 


2.9-8, 


.2- 


-32.4 


30062- 


-02- 


-01- 


-02 


3 


1998/03/26-15:30:25 


712.9-790.2 


1569.2 


32.1±3.8 


27.0 (26) 


2.9-8, 


.2- 


-32.4 


30062- 


-02- 


-01- 


-02 


4 


1998/03/26-17:06:23 


780.6-798.9 


1638.3 


30.6±6.2 


15.7 (12) 


2.9-8, 


,2- 


,32.4 


30062- 


-01- 


-01- 


-00 


1 


1998/03/27-12:27:36 


708.5-743.1 


1332.3 


26.5±4.9 


7.7 (12) 


2.9-8, 


,6- 


,30.9 


30062- 


-01- 


-01- 


-00 


3 


1998/03/27-13:53:23 


723.6-750.5 


1312.7 


34.3±4.3 


15.3 (19) 


2.9-8, 


,6- 


-30.9 


30062- 


-01- 


-01- 


-00 


5 


1998/03/27-15:29:21 


752.0-795.6 


1393.1 


35.2±5.1 


33.8 


(27) 


2.9-8, 


.6- 


-30.9 


30062- 


-01- 


-01- 


-01 


1 


1998/03/28-14:02:52 


769.1-812.0 


1152.2 


28.2±7.3 


22.5 (16) 


2.9-8, 


.6- 


-30.9 


30062- 


-01- 


-01- 


-02 


2 


1998/03/29-09:05:16 


668.1-684.6 


916.6 


36.9±6.6 


5.3 (11) 


2.9-8, 


.6- 


-30.9 


30062- 


-01- 


-01- 


-02 


4 


1998/03/29-10:41:15 


661.5-697.4 


911.2 


39.1±8.6 


17.9 (15) 


2.9-8, 


.6- 


-30.9 
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Table 2. Best fit spectral parameters listed by ObsID and event files, numbered n, within an ObsID. The parameters 
are: the normalization of the disk blackbody temperature kTj„, the normalization Nobb of the disk blackbody 
component, the seed photon temperature kTsee,^, the electron temperature kTg, the optical depth of the 
comptonization cloud r, the line energy of the diskline component, its equivalent width, the 2-20 keV X-ray 
luminosity {L^), the reduced x^, with and without the diskline component. All errors are quoted at the la level. 



ObsID 




n 






^DBB 


^ ' seed 




kTe 




T 




E 




EqW 




Lx 








X2 






X2 




10072 


05 


01 


00 


1 





f.o+0.07 
"°-0.08 


1580.00tiil«,f 





Q7+0.08 
^'-0.08 


2 


66+0 03 
""-0.02 


5 


40+0.13 
^•^-0.11 


7 


271° 


14 
13 


85.46l?2;|5 


9 


01 ± 





01 


29 


54 


(44) 


58 


74 


(46) 


10072 


05 


01 


00 


2 





70+0.07 

'•^-0.07 


1200.28l|?^i? 


1 


nQ+0.08 

"•^-0.08 


2 


66+0-03 
""-0.03 


5 


45+0.13 
^"-0.13 


7 


321° 


18 
17 


76.34li|.«24 


8 


66 ± 





01 


30 


36 


(44) 


53 


11 


(46) 


10072 


05 


01 


00 


3 





71 +0.08 
'1-0.08 


1285.93t|33.07 


1 


09+0.10 
"^-0.10 


2 


6Q+0-04 
"^-0.03 


5 


42+0.15 
^^-0.15 


7 


38i;; 


14 
15 




8 


80 ± 





01 


23 


46 


(44) 


46 


62 


(46) 


30062 


02 


01 


000 


1 





70+0.04 
'"-0.04 


1347.54l«°;f, 


1 


91+0.06 
^1-0.06 


3 


04+0-13 

"^-0.10 


4 


9O+0.26 
^"-0.29 


7 




17 
16 


104.1lli6;i4, 


5 


97 ± 





01 


31 


71 


(37) 


61 


23 


(39) 


30062 


02 


01 


000 


2 





70+0.04 
'■^-0.04 


1 1 c:q c:£:+370.36 


1 


90+0. 06 
^•^-0.06 


2 


Q9 + 0.11 
^^-0.08 


4 


00+0. 24 
'^°-0.28 


6 


991° 


17 
16 


70 Qi +22.86 
10. t>±_i2. 70 


6 


25 ± 





01 


22 


43 


(37) 


38 


50 


(39) 


30062 


02 


01 


000 


3 





71+0.04 
'1-0.03 


1294 09+391-00 


1 


90+0. 06 
^•^-0.05 


3 


Qg+0.13 
""-0.10 


4 


Qg+0.24 
""-0.28 


7 


091° 


18 
18 


76 20+2°-^'* 

(D.^U_jg 21 


5 


93 ± 





01 


27 


60 


(37) 


45 


18 


(39) 


30062 


02 


01 


000 


4 





f;Q + 0.03 

°°-0.03 




1 


99+0.04 
"^''-0.04 


3 


59+0.21 
"^^-0.16 


3 


OC.+0.23 
•^"-0.26 


7 


251° 


14 
13 


103.55ljg g7 


5 


55 ± 





01 


35 


12 


(37) 


67 


42 


(39) 


30062 


02 


01 


00 


1 





f.a+0.03 
°°-0.03 


^A'^9 49+351.69 
l*'J^-^^-268.61 


1 


90+0.03 
^"-0.03 


4 


41+0.29 
^1-0.23 


2 


CO + O.IO 

"■^-0.21 


7 


05lS 


14 
12 


103.2ll?l|° 


5 


35 ± 





01 


42 


08 


(44) 


73 


16 


(46) 


30062 


02 


01 


00 


2 





C7+0.03 
"'-0.03 


1409 04+357.80 


1 


90+0.03 
^°-0.03 


4 


7C: + 0.6S 
'•-'-0.44 


2 


qO+0.31 
•-'"-0.37 


6 


76lS 


14 
11 


127.15l225;05 


4 


24 ± 





01 


57 


42 


(44) 


93 


62 


(46) 


30062 


02 


01 


01 


1 





64+0-03 
"^-0.03 


1757.95l^r,;^I 


1 


90+0.03 
^"-0.03 


5 


O0+0.65 
»U_„ 44 


2 


— O.zo 


7 


061° 


07 
13 


107.87l2f5;™ 


4 


74 ± 





01 


36 


26 


(44) 


70 


81 


(46) 


30062 


02 


01 


01 


2 





(59+0.03 
"^-0.03 




1 


-ifi+O.OS 
l"-0.()4 


5 


60+0-65 
""-0.53 


2 


oc+0.28 
'-'•-'-0.27 


7 


291'^ 


15 
13 


144 1 5+25.85 
l^^-l"-20.59 


4 


71 ± 





01 


43 


39 


(44) 


91 


01 


(46) 


30062 


02 


01 


01 


3 





fil+0.03 
"1-0.03 


2337.68l^»;?« 


1 


17+0.04 
l'-0.03 


7 


9t:+3.9S 
^"-1.05 


1 


7Q+0.34 
'^^-0.74 


7 


43l'> 


19 
19 


115.49li4;2i 


4 


53 ± 





01 


31 


74 


(44) 


56 


95 


(46) 


30062 


02 


01 


01 


4 





fic;+0.03 
""^-0.03 


1641.90t^^r2i 


1 


90+0.03 
^•^-0.03 


5 


Qfi+0.96 


2 


00+0.32 
"•^-0.33 


7 


091° 


16 
14 


103.60l?:^;t| 


4 


62 ± 





01 


50 


44 


(44) 


78 


31 


(46) 


30062 


02 


01 


02 


3 





f.Q+0.03 
"°-0.03 


1264.76tg|;?I 


1 


9C+0.04 
^"-0.03 


3 


45+0.16 

^"-0.13 


3 


07+0.21 
-0.22 


7 


i9ir, 


18 
17 


85.64123:2? 


4 


74 ± 





01 


32 


55 


(44) 


54 


79 


(46) 


30062 


02 


01 


02 


4 





fiS+0-03 
"»-0.03 


1388.04lfo«;?3 


1 


94+0.05 
^^-0.05 


3 


no+0-15 

""-0.11 


4 


OQ+0.29 
"•^-0.32 


7 


25lS 


24 
22 


Q'j Qo-|-24.31 
0' •^•-'-22.82 


4 


94 ± 





01 


59 


54 


(44) 


77 


47 


(46) 


30062 


01 


01 


00 


1 





f.n + 0.03 

"«-0.03 


1142 lQ+283-54 
ll^^-l-^-217.18 


1 


01 +0.03 
■^1-0.03 


4 


fip:+0.92 
""-0.50 


2 


qo+0.37 
•^"-0.48 


7 


221° 


20 
17 


Qfi 42+31-71 
""•^^-18.28 


4 


00 ± 





01 


52 


70 


(44) 


72 


65 


(46) 


30062 


01 


01 


00 


3 





fir:+0.03 
""-0.03 




1 


90+0.04 
^*-0.04 


4 


90+0.59 
^°-0.39 


2 


48+0-35 
^°-0.39 


7 


34lS 


44 
36 


72.18lf,j§ 


3 


93 ± 





01 


42 


92 


(44) 


57 


14 


(46) 


30062 


01 


01 


00 


5 





f;o + 0.03 

"°-0.03 




1 


90+0.04 
^**-0.04 


3 


07+0.18 
•^'-0.14 


3 


qc+0.24 
•^"-0.26 


7 


14lS 


19 
16 


86.96l?°;2^3 


4 


19 ± 





01 


39 


73 


(44) 


61 


13 


(46) 


30062 


01 


01 


01 


1 





fi7+0.03 
"'-0.03 


1048.10l280.i8 


1 


01+0.04 
•^1-0.04 


3 


Of;+0.52 
°"-0.32 


2 


77+0.38 
' '-0.4B 


7 


291° 


11 
11 


133.8012^05 


3 


42 ± 





01 


38 


99 


(44) 


81 


70 


(46) 


30062 


01 


01 


02 


2 





j-o+0.03 
"°-0.03 


1678.761^^9-16 


1 


,7+0.02 
l'-0.02 




5.00 


2 


01+0.02 
•^1-0.02 


7 


571° 


18 
16 


144.39l2f;«i 


2 


68 ± 





00 


51 


30 


(45) 


93 


63 


(47) 


30062 


01 


01 


02 


4 





60+002 
""-0.02 


1471.02t«|;f2 


1 


90+0.02 
^"-0.02 




5.00 


2 


9q+0.02 
^^-0.02 


7 


46lg 


12 
11 


162.04l?3-74 


2 


67 ± 





00 


45 


07 


(45) 


109.00 (47) 
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